Scaling theory of DNA confined in nanochannels and nanoslits 
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A scaling analysis is presented of the statistics of long DNA confined in nanochannels and 
nanoslits. ft is argued that there are several regimes in between the de Gennes and Odijk lim- 
its introduced long ago. The DNA chain folds back on itself giving rise to a global persistence 
length which may be very large owing to entropic deflection. Moreover, there is an orientational 
excluded-volume effect between the DNA segments imposed solely by the nanoconfinement. These 
two effects cause the chain statistics to be intricate leading to nontrivial power laws for the chain 
extension in the intermediate regimes. It is stressed that DNA confinement within nanochannels 
differs from that in nanoslits because the respective orientational excluded-volume effects are not 
the same. 
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A perusal of the rapidly developing literature on 
nanoconfined DNA shows that its behavior is more com- 
plex than anticipated (see e.g. [J 0, 0, H H, @, 0, !, 0, EE 
El El El El EE EE E3 ) • lt appears that more regimes 
are needed besides those originally described by Daoud 
and de Gennes [IH and Odijk . The nanoconfinement 
of a semiflexible chain specifically introduces subtleties 
in the chain statistics which I address here within a scal- 
ing analysis. A complete theory would involve solving 
a Fokker-Planck equation subject to the boundary con- 
ditions arising from nanoconfinement [2fj| . Nevertheless, 
backfolding or hairpin formation may be addressed in a 
mechanical approximation plj though entropic depletion 
of the chain near a wall still needs to be resolved quan- 
titatively [221 ] . Numerical investigations of nanoconfined 
stiff chains interacting via excluded-volume interactions 
have appeared recently [H [24], HE but in the limit of 
ground-state dominance without accounting for hairpin 
formation. 

Let us first consider a very long double-stranded DNA 
chain confined in a nanochannel of square cross section 
whose side D is smaller than the persistence length P so 
that we are in the Odijk regime (D < P). Thus the chain 
may be conveniently viewed as a sequence of deflection 
segments of typical length [Ti| 

A~L>ip3 (I) 

The orientational fluctuations with respect to the channel 
center axis are given by the mean-square average 

^)-Ci(f)" (2) 

It is important to note that the coefficient ci here is quite 
small as has been determined numerically [2?], [28| and 
estimated analytically [12j | . Inevitably, a long chain must 
bear thermally activated hairpins leading to a global per- 
sistence length g as shown in fig I. The entropic depletion 



caused by the nanowalls forces the hairpin bends to be 
tightened up so that g is often considerably larger than 
the persistence length P [2l[. Expressions for g are pre- 
sented in appendix I. 
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FIG. 1: A DNA chain enclosed in a nanochannel of width D. 
The distance between hairpins is typically g. 



In view of the backfolding, segments of the DNA in- 
teract with each other via the excluded-volume effect. 
Owing to the charges borne by the DNA backbone, one 
introduces an effective diameter d e ff rather than a bare 
diameter [2i| (P ^> c? e //)- If the interaction were purely 
isotropic, the excluded volume between a pair of deflec- 
tion segments would scale as [3(| 

P\ - A 2 4// (3) 

But the segments are aligned (see eq J2 J )) so that the 
effective excluded volume becomes [3jJ, [32J, |33| 

E = /3 x (\sm6\) (4) 
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where S is the angle between two deflection rodlets (for a 
computation of the orientational factor, see appendix II). 
It is stressed that we are in the sparse limit d e ff <C D: 
the orientational order is imposed solely by the walls of 
the nanochannel and independent of the density of DNA 
segments. 

The DNA chain of length L may now be viewed as a 
one dimensional walk consisting of L/g statistical seg- 
ments. It is partly self avoiding in the sense that the 
volume exclusion between the deflection rodlets is three- 
dimensional. Employing a mean-field argument of the 
Flory type which is excellent in one dimension [3J], I 
write the free energy of the confined chain as 



Regime 2B: L < L, 

In this case the excluded- volume effect is weak (Z <C 1) 
so the behavior of the DNA is effectively that of an ideal 
chain as long as L ^> g 

Rl ~ Lg (11) 

I now investigate what happens as the nanochannel 
is widened. The global persistence length rapidly ap- 
proaches the usual persistence length (see eq (|25|) in ap- 
pendix I). Accordingly, one may introduce a crossover 
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Here, T is the temperature and ks is Boltzmann's con- 
stant. The first term in eq (JSJ) is the ideal free energy 
needed to extend the DNA chain to root-mean-square 
extension R e and the second term arises from the in- 
teraction of N\ = L/X deflection segments in a volume 
R e D 2 . Minimization of F with respect to R e yields 



= c 2 P 



(12) 



at g ~ P which signals the breakdown of the Odijk regime 
(the numerical coefficient C2 is larger than unity). Con- 
comitantly, the channel is no longer narrow enough to 
impose orientational order on the DNA: ( | sin <5| ) = O (1) 
and E ~ P 2 d e ff. Thus, we now enter the next regime 
upon increasing D. 
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Regime 3: D** > D > D* 
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Eqs ([6]) and ([7]) are reminiscent of those occurring in the 
theory of grafted polymers [35l . 1 3 61 ] . One is now naturally 
led to introduce the following regimes. 



Regime 1: £i > 1 

From eq ([7]) we discern that the chain must be al- 
most fully aligned (R e ~ L) apart from minor fluctua- 
tions given by eq ([2|). 



I now express the total free energy of the chain as in 
eq © but with g = P. This leads to 
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It is again possible to demarcate one subregime in which 
the chain expansion is dominated by the excluded- volume 
effect from another subregime where the chain is ideal 
more or less. The crossover in the contour length is given 
by 



Regime 2A: fi < 1 provided L > L* 

The excluded-volume term in eq ^ is analogous to 
the excluded-volume parameter Z introduced in the two- 
parameter theory of the expansion of flexible polymer 
chains 37] 



N%E_ 
D 2 R e 



(9) 



The excluded-volume effect is fully exerted in the limit 
Z 3> 1. There is a crossover to the case of weakly inter- 
acting segments at Z = O (1). Hence, the contour length 
must be larger than L* if eq |7]) is to remain valid. 
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Superficially, it may appear as if eq (|13p conforms to a 
Daoud-de Gennes type of theory [18| but this is not the 
case for D < D** (see eq (fT6l) ). The intermediate regime 
one has to introduce here is caused by the fact that the 
DNA segments are slender (d e ff <C P). The chain may 
be viewed as a sequence of anisometric blobs, each of 
length (£*P)5 and diameter D. 



Regime 4: D > £>«» 



(10) 



Daoud and de Gennes argued that a flexible poly- 
mer confined in a capillary piles up as a sequence of 



3 



blobs, each blob being viewed as a Flory chain of m 
segments [l^]. The blobs do not interpenetrate owing 
to the excluded-volume repulsion. The supposition is 
that 2> 1 within a blob. In the problem at hand, we 
have m segments of length P interacting by an excluded 
volume P 2 d e ff yielding an excluded-volume parameter 
Z = m 2 f3/m^P 3 = mh (d e ff/P). If the Flory expan- 
sion is to be valid within a blob of radius D, we require 
D ~ [m^ PJ Zi. In other words, upon eliminating m 

5 

we must have Z — (Dd e f f / P 2 ) 6 . Therefore, if the blob 
picture is to be valid, one has to impose Z > 1 implying 
that -D** is expressed by 



£>„ = P 2 /d eff 



(16) 



Note that eq (jTBj) remains valid as can be verified in a 
blob analysis. However, the difference between regimes 3 
and 4 may show up in subtle measurements. 

Let us next turn to nanoslits of rectangular cross sec- 
tion A X D (A > D). In many respects, the reasoning 
is now the same as that presented above so the analysis 
will be brief. I first focus on thin slits (D < irP) in which 
the DNA chain is effectively one dimensional. 

The analogue of eq © is now 
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The orientational factor ( |sin <5|) within E then has a 
rather subtle dependence on A and D (see appendix II). 
The deflection length is still given by eq ([T]). Upon min- 
imizing eq (fl~7|) . we get 
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The second equality in the expression for £3 pertains 
to the limit A ^> D. It is again possible to introduce 
crossovers at £3 ~ 1 and L — L* and so forth. 

If one next increases the width A, the global persis- 
tence length given by eqs (|24|) and (l27|) decreases rapidly 
to the value P at A = C3P where C3 is a numerical con- 
stant larger than unity. At the same time the chain loses 
its local anisotropy: ( | sin <5| ) = 0(1) (see eq (|3ip V It is 
important to realize that the chain remains confined to 
a thin slab (D < ttP). Eq ([TTJ) with g = P still holds 
though £ 3 in eq (fT5| is replaced by 
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Ultimately, if we keep on increasing A, we attain the 
case where the chain may be viewed as a two-dimensional 
pancake (A > R e ). Instead of eq p7|) . we have 
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Minimization of F with respect to R e yields 
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One recognizes the usual 3/4 power law applicable to 
the excluded-volume effect in two dimensions [36]. The 
excluded-volume parameter 
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has to be greater than unity if eq (f2"3"| is to be valid 
otherwise R ~ (LP) 2 . 

The analysis given above has the drawback that the 
numerical coefficients are unknown and may deviate sub- 
stantially from unity. This is exemplified in the appli- 
cation of cqs and (fTU]) to the recent measurements 
on A-phage DNA extended within nanochannels of es- 
sentially square cross-sections Q. For instance, in the 
widest channel of width D = 440 nm, the dimension- 
less parameter £1 is about 0.11 and g — 2.3 /im which 
would lead to a crossover length of about 10 fim (I 
have estimated d e ff to be 4.6 nm on the basis of the 
concentration of buffer used by Reisner et al). But wc 
know that the coefficients in eq J2]) and ([3]) are somewhat 
smaller than unity (see also ref [30(). In addition, there is 
a persistent finite size effect for semiflexible chains which 
significantly suppresses the excluded-volume interaction 
[381 l39j| . On the whole, L* could be an order of magni- 
tude larger so that the DNA of contour length L = 63 
/jm would only be slightly perturbed by excluded volume. 
This would explain why the hairpin theory for the phan- 
tom worm [2l| agrees well with the DNA extensions Q. 
In a similar vein, Krishnan and Petrov [l7| use the same 
theory to explain the ionic-strength dependence of the 
DNA elongation measured in some detail by Reisner et 
al [13j . Theoretically, the strong dependence on salt then 
arises from the exponential dependence of g on the persis- 
tence length P (see eq (1251) ) [17J. The latter quantity was 
dealt with on an empirical level by using the values from 
optical tweezer experiments |4p| . However, Reisner et al 
[13| themselves argue that a blob picture could be valid 
on the basis of assuming g = P (this would be regime 3 
defined above). A reassessment of these experiments is 
warranted to reconcile the apparently opposing points of 
view. 

As I argued above, the degree of orientational order 
depends on whether the cross section of the channel in 



4 



which the DNA is confined is square, rectangular or slit- 
like of infinite extent. In the case of the 1000 x 100 nm 2 
nanoslits used in ref. 12, the quantity £3 = 44 from eq 
(HHJ) turns out to be very large because the ionic strength 
is quite low (d e ff = 79 nm). The DNA is predicted to be 
fully extended even in 100 /xm long slits as is indeed the 
case. Bonthuis at al [l6| have studied the radius of gy- 
ration of DNA in 2d slits as a function of height D right 
into the Odijk regime. There is only one abrupt transi- 
tion at D ~ 2P i.e. there are no intermediate regimes in 
agreement with the analysis presented here. 

In summary, a scaling analysis of nanoconfincd DNA 
has been presented based on orientational order imposed 
by the channel walls and a global persistence length 
greatly enhanced by entropic depiction. Clearly, more ex- 
periments are needed to delineate the regimes proposed 
here. 
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with an optimum radius 
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(E m = 1.5071). In the case of nanoslits, these variables 
are 



F Mt (f) = ^-ln 
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r = (28) 

A + 2E m P v 1 

As A becomes very large, g tends to P, at least to the 

leading order. The limit is not precise because ulti- 
mately fluctuations cause the mechanical approximation 
to break down I21J . 



Appendix II 



Appendix I 

The global persistence length has been computed in 
the mechanical limit for the hairpin configurations (see 
[2lj : the bent DNA remains double-stranded). 



To a good approximation, the orientation-translation 
distribution of the long DNA is that of a chain confined 
in an appropriate harmonic well [12j |. Integrating over 
the translational degrees of freedom, one is left with a 
Gaussian distribution 



g = 3.3082 f exp (F(f)) /k B T 



(24) 



In the case of nanochanncls of square cross section, the 
free energy of a hairpin bend is given by 



F s (f) /k B T 
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for a fluctuating DNA segment (the nanoslit is D wide 
in the x direction and A wide in the y direction; Go — 
{P/D)* and G A ^ [P/A)%, see eq ©). Hence, we have 



2 2 2 2 

(|sin<5|) = J d6 x J d6 y J dw x J dw y f (6 x ,6 y ) f (w x ,w y )\sm5\ 



(30) 



r 



(the appropriate normalization is included in eq (|30p ). where the constant unity has been added to ensure that 
To the leading order integrals like eq (|30[) are readily (|sin5|) = 0(1) in the isotropic limit, 
estimated l3Ul4lj 
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(31) 
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